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I. INTRODUCTION 

In spite of its phenomenological relevance, the infrared 
(IR) limit of QCD has not been fully studied yet because 
of non-perturbative effects that limit the power of stan¬ 
dard tools based on perturbation theory. Our current 
knowledge of the IR limit relies heavily on lattice simu¬ 
lations while, usually, analitycal non-perturbative tech¬ 
niques can only describe the phenomenology by inser¬ 
tion of some free parameters that emerge by some un¬ 
known sector of the theory like vertex functio ns I lH 6fl , 
counterterms PfITTgI| or renormalization schemes [l4|7 A 
mass parameter for the gluon has been shown to cap¬ 
ture most of the non-perturbative effects, leading to a 
reasonable fit of lattice datafl7[-fl9|. but we still miss a 
fully consistent and analytical ab-initio theory without 
spurious fit parameters. 

Quite recently, an optimized perturbation theory has 
been discussed[20|, with zeroth order trial propagators 
that are optimized by some variational ansatz. Many 
variational strategies can be set up by the simple knowl¬ 
edge of self-energy and polarization functions, going from 
the Gaussian effective potential (2llf36| up to Stevenson’s 
minimal sensitivity[37] and the novel method of station¬ 
ary variance [38H40| that has been shown to be viable 
for pure Yang-Mills SU(3) in Feynman gauge (4lj. At 
variance with other analytical approaches, these varia¬ 
tional methods have the merit of reproducing some lattice 
features, like the existence of a dynamical mass for the 
gluon|4l|. without any free parameter, because the trial 
quantities are all optimized by the variational ansatz. Of 
course, the agreement with lattice data is not as good 
as for a fit, but the approximation can be improved or¬ 
der by order and gives an ab-initio description that is 
only based on the original Lagrangian, without spurious 
parameters nor undesired counterterms that would spoil 
the symmetry of the Lagrangian. 

As discussed in Ref. [20], several variational approaches 
can be implemented if the self-energy (the polarization) is 
known, order by order in the optimized perturbation ex¬ 
pansion, as a functional of the trial propagators. On the 
other hand, some internal symmetries of the theory could 
result broken by the truncated expansion and would turn 


out to be only approximately satisfied, so that the actual 
result would depend on further parameters that have to 
do with the gauge choice, the renormalization scheme 
and even the Renormalization Group (RG) invariance. 
All these parameters can also be optimized by a varia¬ 
tional ansatz, yielding an optimal gauge or an optimal 
renormalization scheme [42( |43j. Thus it would be desir¬ 
able to have a general set of explicit expressions for the 
polarization functionals, holding for any gauge, for any 
renormalization scheme and for any trial propagator. Ac¬ 
tually, most of these functionals have been reported for 
a free-particle propagator and in dimensional regulariza¬ 
tion where many terms vanish. A further proliferation 
of terms arises from the use of a generic covariant gauge 
since the trial propagator would be described by two in¬ 
dependent functions for the transversal and longitudinal 
part. Despite of many technical problems, the study of a 
generic linear covariant gauge has attracted some new in¬ 
terest in the last years and the features of the gluon prop¬ 
agator have been investigated on the lattice[Sjj 45| and 
in the framework of Dyson-Schwinger equations jl3l [J0l | . 
Moreover, it has been recently shown that even if some 
IR properties of the gluon propagator, like the dynam¬ 
ical mass, have no effects in the ultraviolet (UV) per¬ 
turbative regime, they can drive a quark-quark interac¬ 
tion that is equal to that extracted by the ground-state 
observables[ 43 , thus enforcing our interest on the gauge 
dependence of the gluon propagator in the IR. 

In this paper we report general integral expressions 
for the ghost self-energy and the gluon polarization, up 
to second order in the optimized perturbation theory, as 
functionals of trial propagators in a generic linear co¬ 
variant gauge, for pure Yang-Mills SU ( N ) theory in any 
space dimension d. The integral expressions hold for 
any renormalization scheme and have been checked by 
a comparison with known results in dimensional regu¬ 
larization and in special gauges like Feynman and Lan¬ 
dau gauge. Then, we use that result for extending to 
Landau gauge a previous calculation of the glu on prop¬ 
agator by the method of stationary variance[20t l38l l39|. 
In fact, the gluo n propagator was studied in Feynman 
gauge in Ref.j41j, while fixed-gauge lattice data are only 
available in Landau gauge. Here, the numerical results 


2 


of the calculation are compared with lattice data in the 
same gauge and with the outcome of the same method 
in Feynman gauge. It turns out that, after renormal¬ 
ization, the gluon propagator is not very sensitive to the 
gauge change and the method seems to be approximately 
gauge invariant, which is a desirable feature of the ap¬ 
proximation. That seems to be in qualitative agreement 
with Ref.[l3]. Some different variational methods are 
discussed and compared, using the same integral expres¬ 
sions for the polarization functionals, but the method of 
stationary variance emerges as the most reliable among 
them. 

The paper is organized as follows: in Section II the gen¬ 
eralized perturbation theory is reviewed and described in 
detail for the special case of pure SU(N) Yang-Mills the¬ 
ory; the first-order graphs for the polarization are evalu¬ 
ated in Section III; the one-loop second-order graphs are 
reported in Section IV, while the two-loop second order 
graphs are evaluated in Section V (the expansion is not 
loop-wise as it is an expansion in powers of the actual 
interaction); In Section VI the gluon propagator is eval¬ 
uated in Landau gauge by the method of stationary vari¬ 
ance and compared with lattice data and with previous 
results in Feynman gauge; a discussion and comparison 
of several variational methods follow in Section VII; some 
details on the numerical integration are reported in the 
appendix. 


A general covariant gauge-fixing term can be written as 


£fix = ^7 Tr 


(d^){d v A v ) 


( 7 ) 


and the quantum effective action r[A'], as a function of 
an external background field A' reads 


e ir [A 1 ] = f D Ae iS[A'+A]j Fp [A'+ A] 
JlPI 


( 8 ) 


where 5[A] is the action, Jpp[A\ is the Faddev-Popov 
determinant and the path integral represents a sum over 
one particle irreducible (1PI) graphsj48|. Since the gauge 
symmetry is not broken and we are mainly interested in 
the propagators, we will limit to the physical vacuum at 
A' = 0, while a more general formalism can be developed 
for a full study of the vertex functk>ns(34j. 

The determinant Jpp can be expressed as a path inte¬ 
gral over ghost fields 

Jfp[A ] = J V^e iS ^ A ’^ (9) 

and the effective action can be written as 

e ir = [ D a ,u,u* e iS ° [j4,u ’“* ] e iSl l A < u > u *] (10) 

Jipi 


II. GENERALIZED PERTURBATION THEORY 

Let us consider pure Yang-Mills SU(N) gauge theory 
without external fermions. The Lagrangian is 


T — Uy m V & fix 

where Lym is the Yang-Mills term 

L Y M = -^Tr (F^F^ 


(1) 


( 2 ) 


and L fi X is a guage fixing term. In terms of the gauge 
fields, the tensor operator F^ v is given by 


where 


F[iv — d^Ay 3 V A^ ig 


A, = Y,X a Al 


A/i, A„ 


and the generators of SU ( N ) satisfy the algebra 


X a ,X b 


= if abc X c 


( 3 ) 


( 4 ) 


( 5 ) 


with the structure constants normalized according to 

fabcfdbc = N5 a d- ( 6 ) 


where the total action in a generic d-dimensional space 
is 

S to t = J £ Y M<l d x+ J £f ix d d x + S gh (11) 

In a generalized perturbation theory[20j, SU we have the 
freedom to split the action in two parts, a trial free action 
Sq and the remaining interaction S'/. We define the free 
action So as 

So = \j A ail {x)D- 1 ab(x,y)A b „{y)d d xd d y 
+ J w„(x)G _1 ab (x,y)u; b {y)d d xd d y (12) 

where Df b v (x,y) and G ab (x,y) are unknown trial matrix 
functions. The interaction is then given by the difference 


Si = S tot - So 


(13) 


and can be formally written as the sum of a two-point 
term and three local terms: the ghost vertex, the three- 
gluon vertex and the four-gluon vertex respectively 


Si — S2 + J d d x [£igh + £3 + £ 4 ] • 


(14) 


In detail, the two-point interaction can be written as 
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$2 — 2 / A afl (x) 


-Do 1 ab(x,y)-D l ab(x,y) A bv (y) d' 


— lM" ( 


where Z} 0 and Go are the standard free-particle propaga¬ 
tors for gluons and ghosts and their Fourier transforms 
read 


' o a b 
Goabip ) = 


p 
Sab 

p 2 


Sab 

T 


D^ h {p) = -^[t^(p) + ^ v {p)] 


(16) 


Here the transverse projector t /IV (p) and the longitudinal 
projector (^(p) are defined as 


Vfp) = 

e r M = Vj f- (17) 

and is the metric tensor. The three local interaction 
terms are 


£3 — ~gfabc{d^A av )A^ > A'^ 

£4 = fabcfadeAb^A cl ,A^A e 

£gh = -gfabdd^DbJbA^. (18) 

The trial functions G ab , Dab cancel in the total action 
Stot which is exact and cannot depend on them. Thus 
this formal decomposition holds for any arbitrary choice 
of the trial functions and the expansion in powers of 
the interaction _Sy provides a generalized perturbation 
theory [2(1 I 40 L HI ], Standard Feynman graphs can be 
drawn for this theory with the trial propagators and 
G ab as free propagators and the vertices that can be read 
from the interaction Si in Eg. ([HI) . As shown in Fig.l, 
we have two-particle vertices for gluons and ghosts that 
arise from the action term S 2 in Ea. lfThl) . while the lo¬ 
cal terms in Eg. (1181) give rise to three- and four-particle 
vertices. The effective action T can be evaluated by per¬ 
turbation theory as a sum of Feynman graphs and sev¬ 
eral variational ansatz can be set up for the best choice 
of the trial functions [20|. mainly relying on stationary 
conditions that can be easily written in terms of self¬ 
energy graphs. Moreover, the propagators can be written 
in terms of proper self-energy and polarization functions 
and their evaluation, up to second order, is the main aim 
of the present paper. First and second order two-point 
graphs are shown in Fig.2. 

Since the propagators are gauge-dependent, we write 
the trial function _D“* as the most general structure that 
is allowed by Lorentz invariance, namely 

D%(p) = Sab [T{p)t^[p) + L{p)^{p)\ (19) 


xd d y + 


u* a {x) [G 0 1 ab (x,y)-G 1 ab (x,y)\uJb(y)d d xd d y 

(Ml 


S 2= + Eq.(15) 

i + y + x Eq <18> 

G - D . - 


Figure 1: The two-point vertices in the interaction S 2 of 
Eq. dlSI) are shown in the first line. The ghost vertex and 
the three- and four-gluon vertices of Eqs. Cl are shown in 
the second line. In the last line the ghost (straight line) and 
gluon (wavy line) trial propagators are displayed. 



(2a) (2b) (2c) 

.-CL- 


(2d) 



Figure 2: First and second order two-point graphs contribut¬ 
ing to the ghost self energy and the gluon polarization. Second 
order terms include non-irreducible graphs. 


while color symmetry ensures that we can always take 

Gab(p) = S ab G(p ) = Sab^- (20) 

p 

where x(p) is a trial ghost dressing function. By the same 
notation, the free-particle propagators in Ea. dlGD follow 
by inserting in Eo. (d) the functions 

t o(p) = —4’ L °(p) = G 0 {p) = -4- (21) 

p2 p2 p2 

Because of the orthogonality properties of the projectors, 
the inverse propagator can be trivially written as 

D-'Zip) = Sab [T(p)-H^(p) + L(jp)-H^(p)\ . (22) 

The trial propagator of Ref.|41] is recovered in Feynmann 
gauge (£ = 1) by taking T(p) = L(p), while in Landau 
gauge (^ —> 0) the longitudinal function L{p) vanishes 
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and the propagator is transverse. In both cases the prop¬ 
agator is described by a single function but in the general 
case two different functions are required. 

III. FIRST ORDER 

Up to first order, the polarization is given by the sum 
of graphs (la) and (16) in Fig.2. The tree-graph II(- la ) is 
just 

- *n (lo) ^ = iD 0 ~Z(p) - iD-'ZiP) (23) 

and in terms of projectors 

n (lo) £» = s ab [uf la) ( P )t^(p) + nf lo) (p)^(p)] (24) 

where 

n fia)b) = t ~\p) +P 2 

n\ la) {p) = L~\p) + P j. (25) 


The one-loop term II(- lt) ) follows from the four-point in¬ 
teraction term £4 in Ea. dlSH that gives the bare vertex 

-plivpa _ - 9 _ rrppvpa . rpppcTV , rppcrvpi 

1 abed L \r~ abed ' ^ aedb ' ^ adbc \ \^J 

where the matrix structure is a product of color 

and Lorentz matrices 

T abcd = RabcdS^ (27) 

with 

Rabcd = feabfecd ( 23 ) 

S p.vpcr = rjPP^vcr _ ( 29 ) 

The one-loop graph (16) then reads 


- * n d b)Z = ^ T aTd I ( 0 j(*^) PWW*) + me^k)] (30) 

and making use of Eq. (JHI we can write 

Ufiftlab = SabNg 2 |(d - l)r? M!/ J ^T(fc) + J ^ [L(k) - T(k)]t^(k) }. (31) 

Integrating in a d-dimensional Euclidean space, for a generic function f(k) that only depends on fc 2 , we can use the 
identity 


/ J 4^j/(M where ^ k ^ = f( k )\w=-k% (32) 


and write the polarization in terms of the constant inte¬ 
grals 

Cn = j ^[T{k E T{k%r 

4 L m = / ^- d [L(k E )T{klr. (33) 

We assume that these diverging integrals are made finite 
by a regulating scheme to be discussed below. The one- 
loop polarization then reads 

n ( i b)Z = -tabV^M 2 (34) 

where the first-order mass term M 2 is defined as 

M2 = Ws^-2) Ko + (d _ 1)/a (35) 


It is useful to introduce the transverse and longitudinal 
massive functions Tm{p), Lm{p) 

PW (p)]' 1 = Po(p)] -1 + M 2 = -P 2 + M 2 

[Lm^p)]- 1 = [LoipT 1 + M 2 = -^- + M 2 (36) 

and the massive propagator 

D m ^(p) = T M (p)t^(p) + L M (p)e V {p) (37) 

that describes a free massive particle in a generic co¬ 
variant gauge. In the special cases of Feynman gauge 
(£ = 1) and Landau gauge (£ —> 0) we recover the mas¬ 
sive free-particle propagators D^ip) = V^/(—p 2 + M 2 ) 
and DZ(p) = t^(p)/{~P 2 + M 2 ), respectively. With 
that notation, the total first-order polarization III can 
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be written in the very simple shape 

n io6 = n (l a)Z + n (16)ab = ^ ab ~ <^M ^ (38) 

There is just one first-order graph for the ghost self¬ 
energy, arising from the two-point non-local term in 
Ea. (115D as shown in Fig.2, so that the first-order self¬ 
energy can be written as 


According, in Feynman gauge the larger effective cou¬ 
pling would give a larger mass as compared with Landau 
gauge. 


T, a 1 \p) = 6 ab [G- 1 (p)-Gc\p)\. (39) 

The Gaussian Effective Potential (GEP)[HH3i| can be 
derived by the requirement that the functional derivative 
of the first-order effective potential with respect to the 
trial functions D and G is zero, that is equivalent[20] to 
the self-consistency condition of a vanishing first-order 
self-energy and polarization, III = 0 and Ei = 0. The 
gap equation that arises has been first investigated by 
Cornwall [41] in 1982 as a simple way to predict a gluon 
mass. In the present formalism, from Eas. (l38l) and (l39l) 
the stationary conditions for the GEP, that derive from 
the vanishing of first-order self-energy and polarization, 
give a decoupled ghost with G = Gq and a free massive 
gluon with D = Dm , where the mass M follows from the 
gauge-dependent gap equation (13511 that can be formally 
written, by a change of argument in the second integral, 


N(d — l) 2 ff| r d d k E 1 
d J (2n) d k 2 E + M 2 


(40) 


where the gauge dependence has been absorbed by the 
effective coupling 


2 2 
9$= 9 


1 + 


(d /2 ' 

d- 1 


(41) 


IV. SECOND ORDER - ONE LOOP 


The generalized perturbation theory that arises from 
the expansion in powers of the interaction Si is not a 
loopwise expansion in powers of the coupling constant, 
so that one-loop and two-loop graphs coexist in the sec¬ 
ond order term of the polarization. The standard one- 
loop graphs, namely the ghost and gluon one-loop graphs, 
II( 2 0 ) and II( 2 b) in Fig-2, are described in this section to¬ 
gether with the one-loop ghost self-energy. The other 
second-order terms, namely the second-order one-loop 
graph II( 2 (Z) and the two-loop graphs II( 2 C ) and II( 2 e ) will 
be discussed in the next section. 


A. One-loop ghost self-energy 


The one-loop ghost self-energy follows from the bare 
vertex of the ghost-gluon interaction term L g h in Ea. lfTSl) 


tE a ^(p) — g fcdafc'bd' 


d d k 


kft)p v [iD^,(k)iG dd '{p - k)\. 


(42) 


Making use of Eq.([G]) and integrating in a d-dimensional Euclidean space, we can split the self-energy in two terms 


S 0 b(p) = Sab [£ T (p) + £ L (p)] 


(43) 


where 


£ t (p) 

£ l (p) 



xi^E ~PE)T(k E ) 
(k E — Pe) 2 

x( k E -PE)L(k E ) 

{k E - p E ) 2 


2 ( k E -PE ) 2 ’ 

Pe p 

^E J 

( k E -PE ) 2 n 4 

-T2- i k E ■ Pe) 


(44) 


These integrals are functionals of the ghost dressing function \ of Eq. d^Ul) and of the gauge-dependent transverse and 
longitudinal gluon propagators, respectively. 
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B. Ghost loop 

The one-loop polarization term II( 2 a)> the ghost loop in Fig.2, also follows from the bare vertex of the ghost-gluon 
interaction term L gb in Eg. (1181) 


- *n (2o) £(p) = ~g 2 f ab chad J p + kYk v iG{p + k)iG(k). (45) 

Making use of Eq. (© and integrating in a d-dimensional Euclidean space, with the same notation of Eq. (|24p , we can 
write the transverse and longitudinal parts in terms of the trial ghost dressing function % of Eq. (1201) 


n 02a) (?) 

n (2a) (P) 


Ng 2 

~W~T) 

- N »i 


f d d k E x(pe + k E )x(k E ) 

1 

(k E ■ p E ) 2 

J ( 2 n) d (p E + k E ) 2 


k\v\ . 


d d k E x(pe + k E )x(k E ) 
( 27 r) d ( p E + k E ) 2 k 2 E 


(k _e • Pe) + 


(k E ■ Pe) 
Pe 


21 


in agreement with the result reported by other authors 


0. 


(46) 


C. Gluon loop 


The one-loop polarization term 11(25), the gluon loop in Fig.2, follows from the gluon-gluon interaction term L 3 in 
Ea. (ll 8 D that gives the bare three-particle vertex 


Kbc 0 > k) = ^ {YYP P - <f) + vTifY - n + rf<>{k v - p v )} . 


(47) 


The one-loop graph (2b) then reads 

_, ;TT fn/ f \ _ 313! /' d d k 

1 ( 2 b ) oa'W 2 J ( 2 n) a 


■iD bb ',pp'(jp + k)iD cc ^ TT f(k)T^(p, -p - k, k)Y v J b J c ,(-p,p + k, -k) (48) 


and since the trial propagator is defined by two indepen¬ 
dent functions, there are 36 terms for each of the lon¬ 
gitudinal and transverse parts of the polarization. We 
can write them in a more compact shape by introducing 
some degree of redundancy in the notation. By Eq.©, 
the sum over color indices gives a diagonal matrix, so 
that we can use the same notation of Ea. (l24l) and drop 
all color indices. Let us denote by a, /?, 7 the three mo¬ 
menta in the vertex, a = —p, (3 = p + k, 7 = — fc, so that 
a + (3 + 7 = 0 . Then, denote by , Bg v , C£ v the three 
projectors 

A^=PoT{a) 

B^ v = PY v (p) 

CY = P/"( 7 ) (49) 

where a, b, c = ± 1 , while P± are the transverse and 
longitudinal projectors P£ v (k) = t^ v (k) and P-^ v (k) = 
i^ v (k), that can also be written as 

PY ( k ) = n a Y v - ai» v (k) (50) 

where n a = (1 + a)/ 2 . 

Moreover, with the same notation of Ea. (H?41) . let us 
denote by A ai 23b, C c the numbers 

A + = (d - l)nf 26) (a), A- = n^) (a) 

23+=T(/3), 23_ = L(/3) 

e+ = r( 7 ), e_ = l( 7 ) 


so that, having dropped color indices (not to be confused 
with the sign indices a, &, c in this section), we can write 


A a = AYn ( 2 % j P ) 

D^( P +k) = 

6 =± 1 

D^(k) = e cCY- (52) 

c=± 1 


Inserting in Ea. (l48l) . the transverse and longitudinal po¬ 
larizations A a can be simply written as 


A 


a 



id d k 

(2Tr) d 


3 abc (a,/3, 7 ) 


( 53 ) 


where, with the obvious shorthand notation 


(51) 
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k ■ X • ■ • Y ■ p 


= k li X»---Ylp v , 


X ■ Y 


= X tJV Y, 


Vf! 1 


(54) 


the matrix IF reads 

J a bc(o:,l3, 7 ) = (a — (3) ■ C c ■ (a — /3) 


An ■ Bh 


+ 2 


(a - /3) • C c ■ B b ■ A a ■ (/3 - 7 ) + cycl. perm., 


(55) 


summed over the three cyclic simultaneous permutations of all the arguments, indices and projectors, i.e. a —> (3 
7 — > a together with a —> b —> c —> a and A—> B —> C —>■ A. A straightforward but tedious calculation yields 


?abc(a,/3, 7 ) = 3 (n a n b n c )(d- 1 ){a 2 + /3 2 +'y 2 ) + { ( acn b )/3 


1 - 


(/3 2 — a 2 — 7 2 ) 2 


4a 2 7 2 


cycl. perm. 


( an b n c ) 


a 


-(/3 2 + 7 2 )-(2d-3) 


Q3 2 ~ 7 2 ) 2 
2 a 2 


+ cycl. perm. 


(56) 


The result holds for any gauge parmeter £, space di¬ 
mension d , and for any trial propagator. It has been 
found in a gree ment with previous calculations in Feyn¬ 
man gauge [44| and with older results for a free-particle 
trial propagator in generic covariant gauge[50], 

For instance, in Feynman gauge, the polarization func¬ 
tion II(p) of Ref. (41j is defined as II(p) = {A+ + AJ)jd 
and the trial propagator is taken to be T{k) = L(k) so 
that the kernel of the integral in Ea. (153l) is just given by 


^ ^abda, 13, 7 ) = 6( ' d d 1 \ p 2 + k 2 +p- k ) (57) 

abc 


in agreement with Eq.(All) of Ref. fdA~l| for d = 4. In the 
same work the function lF(p) is the transverse polariza¬ 
tion n'(p) = A+/(d — 1) and the corresponding kernel in 
Eg. (1551) . for d = 4, is given by 


—y PWa, 13, 7)L=i = 5 P 2 + 2 fc 2 + 2p-k 



(p-k) 2 ' 
k 2 p 2 


(58) 


in agreement with Eq.(Al2) of Ref.jSJ. In the work 
of Watson[5(il| the trial propagator is taken to be the 
gauge-dependent free-particle function Dq as defined in 
Eas. lfTfil) . (12T1) . The function Jp 1 ' in Eq.(3.3.4) of that 
work[H(| is defined as (d — 1 (Jp 1 ^ = (A + +A-)/p 2 and 
the corresponding kernel in Eg. (1531) can be written as a 
polynomial [wq + + W 2 ^ 2 ] with coefficients 



a 


[9 r afcc(o!, f3, 7 )]6= c =l 



^2l[%bc(a,P, 7 )] 

a ^ 


6 =+i T [T a b c (ck, j3, 7 )] 

c— — 1 


b=-l 

c=+l 



a 


[A a bc(. a , f3, 7)]b =c= — 1 ■ 


(59) 


The coefficients can be easily evaluated by Eq. (1561) and if 
we drop all terms that vanish under integration (because 
of symmetry or by dimensional regularization) we obtain 


7 a 2 a 

W0 =3d-- + 2(2-d) 7 - i; ^ 


of 


W\ = 1 + (2d - 5)— + 
T 


a 


2j 2 f3 2 


1 a 2 a 4 
= 2 + t 2 ~ 4j 2 /3 2 


(60) 


in agreement with Eq.(3.3.4) of Ref.[50| that was evalu¬ 
ated by a computer routine for algebraic computations. 
The general result in Eq. (157)1) holds for any choice of the 
trial propagator and contains all terms that might not 
vanish by symmetry or dimensional regularization when 
a generic massive propagator is considered. Moreover the 
result does not depend on a specific regularization scheme 
and can be used for any kind of calculation. 

In Landau gauge, the propagator is transverse and is 
defined by only one function T (p ), so that the transverse 
polarization A+ is obtained by retaining only one term, 
for b = c = +1, in Eo. H53D . and the corresponding kernel 
for d = 4 reads 


[Aabc(o!,/3,y)\ a _ b _ c _ +1 — 1 

x ll(fc 2 + p 2 ) + 2(k ■ p) + — 


( k-p ) 21 

k 2 p 2 

4 10 p 2 k 2 + k 4 ) 


{k+p) 2 


■ ( 61 ) 


V. SECOND ORDER - TWO LOOP 


Besides the standard one-loop graphs of the previous 
section, the second-order polarization includes the one- 
loop and two-loop tadpoles n ( 2 ( 2 ), Il^e) an< l the two-loop 
sunset II( 2 c)- 

































A. Tadpoles 

The one-loop and two-loop graphs, H( 2 d) and II( 2 e) in 
Fig.2, follow from the standard tadpole II( lfc ) by insertion 
of the total first-order polarization in the loop 


The constant integrals /J m . I% m were defined in Eg. ( 1551 ) 
and are functionals of the transverse and longitudinal 
trial functions T(p), L(p ), respectively. 


iD^ iD^{-iIh pa )iD av (62) 

that is, making use of Eq. (l38l) , 

D' lv D^ D av . (63) 

Insertion in Eas. (f34ll . (HTSll yields 


B. Two-loop sunset 


n (2d) ^ + n (2e) ^ = -«- [m 2 - m 2 ] 

where the new mass constant reads 


M 2 = N ^d-l) 
d 


T l 

i 2,l 

e 


(64) 


o + ^ + (d- 1 )(M% 0 + 1^) 


(65) 


The two-loop graph (2c) in Fig.2 is the most involved 
and even if its contribution is small when the coupling 
is not too large, it can be relevant in a variational ap¬ 
proach when the coupling is allowed to increase enough. 
The calculation follows from the explicit four-gluon bare 
vertex in Eo. (f26ll 


J 


-in 


(2c) 


MM 


4!4! 


1 f'n'l — vM^pa- t^m v p a 

' tel — 0 | 1 abce 1 a'b'c'e' 


d d k 

(27r) d 


d d q 

(2n) d 


d d t 

(27r) d 


[iDbb>,uv'{k)\ [iD cc t' PP >(q)\ [iD ee i^ a >{t)\ (2n) d S d (k+q+t+p) 


( 66 ) 

and the explicit symmetry by permutation of dummy integration variables ensures that, under integration, the whole 
expression is symmetric for the exchange of the corresponding Lorentz and color indices in the matrix factors. Thus 
the sum over permutation in Ea. (|26ll can be replaced by a factor of three in one of the vertices, yielding 


n 


(2c) 


MM' ( n ) _ 9_rpnvpa 
nn' ) o ^ abce 


rriH V p G . rpP p G V . rpP G V p 
n'bce ^ n r r&h ' 


ceb 


' ebc 


E 

ijk=±l * 


id d k 


id d q 


(2n) d J 


B° 




(67) 


with a compact notation that extends that of the previous 
section: here we define the four vectors a = p, /3 = 
—(p + q + k), "f = q, e = k so that a + /3 + 7 + e = 0 and 
add a new projector to the set in Eg. ( 147 ) 1 ) 


E^ = P e > iv {e). ( 68 ) 


Moreover, in this section, we denote by A ai CB&, C c , £ e the 
numbers 


A+ = (d- l)nf 2 c) (a), A _ = nf 2 c) (a) 

® +=T(p), V 

e+ = T( 7 ), e_ = L( 7 ) 

£+ = T(e), £_ = L(t) (69) 


so that, dropping color indices, we can write 
Aa = ir n (2c)fti » 

i>'“'(p+9+fc)= 53 

6=±1 

n^(g) = ^3 e c cr 

c=± 1 

D^ik) = 53 £c£r- (70) 

e=±l 

Under integration, the matrix structure in Eg . (1671) sim¬ 
plifies because of the permutation symmetry of dummy 
integration variables and the three matrix products can 
be recast as a single Lorentz matrix that multiplies 
three color matrices 

H(2c)((q/ = [2 Ra'bceRa'bce Ra'bceiRa'ceb 4” Ra'ebc)\ 

= 3N 2 S aa (71) 

where the last equality folows by Jacobi identity. Then, 
dropping color indices, the transverse and longitudinal 
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polarizations A a can be written as 


A 


a 


3 N 2 g i 
2 



id d q 

(2ir) d 


® & e c £ e J obce (a,/3,7, 


e) 


(72) 

where the kernel T follows from the projection of 
by the projector A according to Eg. (171)1) and with the 
shorthand notation of Eg . (1541) it can be written in terms 
of traces of projectors 


Because of the obvious symmetry of the integral in 
Eo. (l72l) . the result is invariant for any simultaneous per¬ 
mutations of the last three arguments, indices and projec¬ 
tors, i.e. /? — > 7 —>■ e p together with b —> c — > e -» b 
and B —»• C —> E —>• B. Using that symmetry the kernel 
If can be written as 


■Eabce (ct, , 7; c) 


A a ' E e 

Cc-B b 

— 

A a • E b ■ C a ■ B b 




(73) 


J 


Aabce{a, P, 7, e) = d(d - l)(n 0 n b n c n e ) - 3 (d - 1 ){n a bn c n e ) - (d - l)(an b n c n e ) + ( n a n b ce ) 


(e * y ) 2 

2 +( d -3rS- 


( an b cn e ) 


(d~ 3) 


(a- 7 ) 


21 


a 2 7 2 


- ( n a bce ) 


(e-7) 

e 2 7 2 


(e-7) - 


(P ' l)(P' e ) 

A 2 


(an b ce ) ^ ^ 


6 7 


(e' 7) — 


( 7 
(a • 7 ) (a • e) 


2 (abn c e) 


(a- P) 

a 2 p 2 


{a-/3)~ 


(a-e)(/3-e) 


(abce) ( a f [(a • P)("f • e) - (a • e)(/3 • 7 )]. 

I- 


(74) 


That explicit expression gives the two-loop sunset graph 
(2c) by a double integration in Eo. flTTl) and holds for any 
covariant gauge, any trial propagator and any space di¬ 
mension. For instance, in Feynman gauge, the polar¬ 
ization function II (p) of Ref.[44| is defined as II (p) = 
(A+ + A-)/d and the trial propagator is taken to be 
T(k) = L(k) so that the kernel of the integral in Eo. flTTl) 
is a constant and is just given by 

^J2^»ce(a,p, 1 ,e) = (d-1) (75) 

abce 

in agreement with Eq.(Al9) of Ref.j4l| for d = 4. In the 
same work the function n'(p) is the transverse polariza¬ 


tion II ' (p) = A+/(d — 1) and the corresponding kernel in 
1x1. (1771) is given by 

J—jYl [ J ^e(a,/3,7,£)] a= i = (d- 1) (76) 

bee 

again in agreement with Eq.(Al9) of Ref. 0. 

In Landau gauge the result is more involved. The prop¬ 
agator is transverse and is defined by one function T(p), 
so that the transverse polarization A+ is obtained by re¬ 
taining only one term for 6 = c = e = +lin Eo. (l72D . 
The corresponding kernel for d = 4 reads 


[Aabce(a, P, 7, e)] Q= fe =c= e=+l = 4 + 
( k-q ) 2 


(k-q) 


+ 


a=fe=c=e=+l * 1 jAq 2 

( P-q ) 2 ,2 [p- (p + q + k )] 


[q ■ (k + q+p)} [k- (k + q + p )] (p ■ q)(p ■ k) 


(k + q + p) 2 


k 2 q 2 p 2 q 2 p 2 {p + q 
(k-q)[p-(p + q + k)} 


k) 2 


p- (p + q + k) — 


(p ■ k) [k ■ (k + q+ p)] 


k 2 


p 2 q 2 k 2 (p + q + k ) 2 


[p-(p+q + k)\ (k -q)-(p-k)[q-(k + q + p)\ 

I- 


(77) 


VI. STATIONARY VARIANCE IN LANDAU 
GAUGE 

Explicit expressions for the second-order graphs are 
useful for a direct comparison of variational results in 


different gauges. For instance, the method of station¬ 
ary variance [20|, [38H40|J has been shown to be viable in 
Feynman gauge where provides a good description of the 
gluon propagator[4lj. Here, we explore the outcome of 
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the same method in Landau gauge where lattice simula¬ 
tions are available [Hll. f52| . 

The variance is stationary when_the trial propagators 
satisfy the stationary conditions|20] 

n 2 = iii 

E 2 = Ei (78) 


where Ei is the first-order ghost self-energy in Ea. (l39l) . 
E 2 is the sum of reducible and irreducible second-order 
graphs, Iii is the sum of the first-order polarization 
graphs (la), (16) and II 2 is the sum of reducible and 
irreducible second-order polarization graphs. With the 
notation of Eas. (l38l) . (|39ll . (l43l) the stationary equations 
read 


J 


n—i 

^ ab 


- 8 ab D- M l ^ = [D - 1 Z - S ac D - 1 w] [D C e, p<t] [ D- 1 % - 8 eb D^ + II* £ 

G~ l - Gq 1 = [G- 1 - G^ 1 ] G [G- 1 - G^ 1 ] + E T + E L 


where the proper polarization Lf 2 is the sum of all the irreducible second-order terms 

TT* M*' _ TT | TT . T | fllS . tt 

-^^(2a) ab ^(2 b) a b ^^(2c) a b ^(2 d) a b 


± 2 ab 


nr** Z + n 


fiu 
(2e) ab ' 


r 


(79) 


(80) 


The coupled set of integral equations (1791) can be written 
as 


D Z = 8 abD M - Dm w n* ab , P oD M ° v 

G = Go - Go [£ T + E L ] Go (81) 

and hold for any gauge. The first of Eqs. m shows that 
the optimal propagator D Z must be diagonal in color 
indices. 

In Landau gauge (£ —> 0), according to Eo. (l37ll the 
massive propagator Dm M1/ becomes transversal so that 
any longitudinal term in the polarization does not play 
any role in the first of Eas. dSlTl . yielding a pure transver¬ 
sal solution for the optimal propagator D Z- Moreover 
in the second of Eqs. (15TT) the longitudinal term £ L is 
zero according to Eq. (J44]) and we can drop it, yielding 
a decoupled set of integral equations for the transver¬ 
sal components. With the same notation of Eg. (1241) the 
transversal component of the proper polarization II£ in 
Eg. (1801) can be written as 

nf Z(p) = ^mf(p) (82) 

where the scalar function (p) is the sum of all the 
transversal second-order irreducible polarization graphs 
(2a), (2b), (2c), (2d) and (2e). The optimal propagator 
can be written as 


DZ(p) = S ab t^(p)n P ) (83) 

where the scalar function T(p) satisfies with G(p) the 
coupled set of stationary equations (15TT) that in Landau 
gauge become 


T{jpe) = ^5 

G(p E ) = 


1 


Pe + M 2 
1 


1 - 


n ?{pe) 


Pe L 


1 - 


Pi + M 2 
X t (pe) 

-P 2 e 


(84) 


and we have made use of the explicit expressions of Go 
and Dm in the Euclidean formalism. The last equation 
can be written in terms of the dressing function Eg. (1201) 
as simply as 

x(pe) = 1 + (85) 

Pe 


The gauge-dependent mass parameter M is given by 
the one-loop graph (lb) and follows from its definition 
in Eo. (l35l) which closes the set of equations and must 
be evaluated by insertion of the actual propagator T(jp) 
instead of the first-order massive propagator that was 
used in Ea. (l40l) for the GEP. Of course, the mass pa¬ 
rameter does depend on the choice of gauge, as it was 
obvious at first order. Thus, it defines an energy scale 
that must be not confused with the actual gluon mass of 
the renormalized propagator. Some typical values of the 
mass parameter are reported in Table I and compared 
with the corresponding values in Feynman gauge as dis¬ 
cussed below. In Landau gauge the gap equation (1351) 
reads 


M 2 


Ng 2 (d — l) 2 tT 


( 86 ) 


where the integral 7^ 0 is defined in Eq. (1331) and is a func¬ 
tional of the unknown full propagator T(p). The station¬ 
ary equations (1571) together with the gap equation (157)1) 
can be made finite by a proper regularization scheme and 
solved numerically. Details on the numerical calculation 
for d = 4 are reported in the Appendix. 

In this work the integrals are regularized by a finite 
cutoff A in the Euclidean four-dimensional space (d = 4) 
where we take p 2 E < A 2 . The simple choice of a cutoff 
gives physical results that are directly comparable with 
the outcome of lattice simulations where a natural cut¬ 
off is provided by the lattice spacing. The bare coupling 
g = g (A) is supposed to depend on the energy scale A and 
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g M l (MeV) M F (MeV) 
0.5 172 457 

1.0 187 573 


Table I: Mass parameter in Landau gauge (Ml) and in Feyn¬ 
man gauge ( Mf ) for d = 4 and N = 3. 



P 2/ A 2 g=! 


Figure 3: Log-log plot of the renormalized propagator Tr(p) 
as obtained by appropriate scaling of the bare propaga¬ 
tor for N = 3, d = 4 and for a bare coupling g = 
0.50,0.60,1.00,1-20,2.80. The scale is arbitrary because of 
scaling: all curves have been scaled in order to fall on top 
of the g = 1 bare propagator. Energy is in units of A 3= i so 
that g( 1) = 1 (for g = 1 the curve is not rescaled). Some de¬ 
viations from scaling become more evident at the very large 
coupling g = 2.8. 

RG invariance requires that the physical content of the 
theory should be invariant for a change of scale A —»• A' 
followed by a change of coupling g —> g'. Then, physical 
renormalized functions that do not depend on the cutoff 
can be obtained by scaling. It is important to point out 
that the present regularization scheme does not need the 
inclusion of any counterterm in the Lagrangian and espe¬ 
cially mass counterterms that are forbidden by the gauge 
invariance of the Lagrangian. The interaction strength g 
at a given scale A is the only free parameter while the 
function g( A) can be determined by RG invariance. In 
principle, one could fix the scale by a comparison with 
experimental data. However, in the present model cal¬ 
culation on pure Yang-Mills theory, we will fix the scale 
by a comparison with the lattice data. Since the origi¬ 
nal Lagrangian does not contain any scale, it is useful to 
take A = 1 and work in units of the cutoff, at a given 
bare interaction strength g. Thus the choice of A will 
be equivalent to a choice of the energy units. RG invari¬ 
ance requires that a renormalized propagator Tr(p) can 
be defined at an arbitrary scale p by multiplicative renor¬ 
malization, that is equivalent to say that by scaling all 
bare functions at different couplings can be put one on 
top of the other. Of course, since the approximations and 



p 2 (GeV 2 ) 


Figure 4: The renormalized propagator T(p) in Landau gauge 
for IV = 3, d = 4 and for a bare coupling g = 0.5. The 
scale has been fixed in order to fit the lattice data of Ref. [51] 
(g = 1.02, L=96) that are displayed as filled circles. The 
propagator in Feynman gauge is shown for comparison as a 
dotted line. 



Figure 5: The renormalized propagator T(p) in Landau gauge 
for N = 3, d = 4 and for a bare coupling g = 1.0. The 
scale has been fixed in order to fit the lattice data of Ref. na 
(g = 1.02, L=96) that are displayed as filled circles. The 
propagator in Feynman gauge is shown for comparison as a 
dotted line. 

the numerical integration spoil the scaling properties, we 
will consider the scaling as a test for the accuracy of the 
whole procedure. We can study the scaling behaviour in 
a log-log plot where the bare functions should go one on 
top of the other by a simple translation of the axes. 

In Fig.3 the renormalized gluon propagator is shown 
for N = 3, d = 4 and for several couplings in the range 
g = 0.5—2.8. Scaling is rather good in the range g = 0.5— 
1.2, but gets spoiled at the rather large coupling g = 2.8. 
That could be a limit of the second-order approximation. 

For any coupling, the energy scale can be fixed by com¬ 
parison with the lattice data, yielding a physical renor- 
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Figure 6: The total polarization is displayed together with the 
second-order contributions of the one-loop graphs (2a), (2b) 
and of the two-loop graph (2c) for N = 3, d = 4. The coupling 
is g = 1 (/3 = 6). The total polarization includes the constant 
contributions II( 2 d) = —5.25 ■ 1CP 3 and II( 2 e ) = 2.87 • 10 -3 in 
units of the cutoff. 



p 2 /A 2 


Figure 7: Same as Fig.6 but for a strong coupling g = 2.8. The 
constant terms are II( 2 d) = — 7.30-10 -2 and II( 2 e) = 3.48-10 -2 
in units of the cutoff. 

malized propagator T{p) that is shown in Fig.4 and Fig.5 
for g = 0.5 and g = 1.0 respectively. 

As shown in Table I, when expressed in physical units 
the mass parameter M is almost constant with respect to 
changes of the coupling, while it remains very sensitive to 
the choice of gauge. However, after scaling, we can define 
a physical mass m 2 = T(0) _1 that does not depend on 
scale and gauge because it is made to coincide with the 
lattice value m ~ 320 MeV. 

A direct comparison of results in Feynman gauge (41) 
and Landau gauge shows that the renormalized propaga¬ 
tor is not very sensitive to the choice of gauge, confirm¬ 
ing that gauge invariance survives albeit approximately. 
Thus the effects of a different mass scale, more than dou¬ 
ble in Feynman gauge, are absorbed by renormalization. 


As shown in Fig.4 and Fig.5, where the propagator in 
Feynman gauge is reported for comparison, we can say 
that the difference between the lattice data and the cal¬ 
culated curves cannot be ascribed to gauge differences 
but is probably a consequence of the finite order (second 
order) of the approximation or of the finite size of the lat¬ 
tice. In fact, even if optimized by a variational method, 
the nature of the calculation is perturbative and can be 
improved by inclusion of higher orders. Actually, we do 
not expect that a perfect agreement could be reached in 
the UV limit of standard perturbation theory because of 
the simple renormalization scheme of the present calcu¬ 
lation that is based on an energy cutoff. In that scheme, 
the spurious quadratic divergence would spoil the result 
for the UV limit. The problem of cancelling that diver¬ 
gence without affecting the IR limit has been discussed 
by several authors and recently reviewed in Ref. 0 . It 
is a major problem that has not found a satisfactory so¬ 
lution yet. 

Finally, it is instructive to compare the weight of the 
single graphs in the second order polarization for the 
studied case N = 3, d = 4. At a rather strong cou¬ 
pling g = 1 (/3 = 6), the total polarization is reported 
in Fig.6 together with the contributions of the ono-loop 
graphs (2a), (2b) and of the two-loop graph (2c). The 
total polarization includes the constant contributions, 
II(2d) = —5.25-10 -3 and n( 2e ) = 2.87-10 -3 in units of the 
cutoff, that are not negligible compared to the one-loop 
graphs. We observe that the two-loop term (2c) is very 
small and rather constant so that it could be neglected 
without serious consequences. On the other hand, at 
the very strong coupling g = 2.8, Fig.7 shows that the 
two-loop graph is still rather constant but becomes quite 
important in the low energy limit where it is as large 
as the total polarization. The constant terms are also 
rather large and amount to n ( 2 d) = —7.30 • 10~ 2 and 
II(2e) = 3-48 • 10 -2 in units of the cutoff. 


VII. DISCUSSION 

By the explicit knowledge of the second-order polariza¬ 
tion, several variational strategies can be set up for the 
optimization of the perturbation expansion. The method 
of stationary variance [2^, l38l - l4l| has been discussed in 
the previous section. The generalized perturbation the¬ 
ory can also be optimized by other methods like minimal 
sensitivity[373 or by the self-consistent requirement of a 
vanishing self energy. Here, we give a brief description 
and comparison of some different methods. 

As discussed in a recent paper (20], the stationary con¬ 
ditions for Stevenson’s method of minimal sensitivityj37] 
can be written as simply as 

n 2 = 0, E 2 = 0 (87) 

to be compared with Eas. flTSl) for the stationary vari¬ 
ance. Eas. (l87l) are equivalent to the requirement that the 
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Figure 8: The self-consistent renormalized propagator T(p) 
in Landau gauge for IV = 3, d = 4. The scale has been fixed 
in order to give a rough fit of the lattice data of Ref. [51] (g = 
1.02, L=96) that are dis play ed as filled circles. The second- 
order propagator of Ref.[4jj in Feynman gauge is shown for 
comparison as a dotted line. 


second-order effective potential is stationary for any vari¬ 
ation of the trial propagators. We explored this method 
in Landau gauge, but Eas. (187l) have no physical solution 
for N = 3 and d = 4. In fact it is well known that 
sometimes that method does not show any range of pa¬ 
rameters where the effective potential is stationary. A 
second derivative would be required for imposing that 
the sensitivity is minimal. 

An other simple approach would be based on a 
first-order optimization of the expansion, followed 
by a second-order evaluation of polarization and 
propagator[53]. That would be equivalent to taking 
the trial propagator equal to the massive propagator 
Dm = (—p 2 + M 2 ) -1 . Even if the mass would not be 
dynamical in the trial propagator, a second order propa¬ 
gator D 2 can be defined as usual by Dyson equations 

D 2 - 1 (p) = D- 1 (p)-n 1 (p)-n*(p) (88) 


where D(p) is the trial propagator and D = Dm in the 
actual case. Since the first-order optimization is self- 
consistent, it requires that II1 = 0 and in Landau gauge 
by Eos. (I82p . (1831) we can define a transversal second-order 
propagator T 2 that reads 


t 2 (p) 



(89) 


Thus the second order propagator would acquire a dy¬ 
namical mass that in the low energy limit tends to 
m 2 = M 2 — n 2 T (0). The advantage of this basic approach 
is that the trial propagator is simple and the theory can 
be renormalized by standard dimensional regularization. 
A similar massive model has been recently studied p~7l— 
El and shown to be in close agreement with the lattice 
data. However the mass was regarded as a free param¬ 
eter rather than a variational parameter. It would be 


interesting to see how close the result would be for the 
variational approach that does not contain free parame¬ 
ters at all. 

The simple first-order optimization would not be self- 
consistent at second order but one can attempt and ex¬ 
tend it by a self-consistent approach. Ea. (l88D is quite 
general and can be made self-consistent by the simple 
requirement that the total proper polarization vanishes 
exactly 


n 1 ( P ) + n$(p) = o. (90) 


That would generalize the first-order stationary condi¬ 
tion ni = 0 which holds for the GEP . If the polar¬ 
ization were not truncated at the second order, Eq. flUUl) 
would be the exact condition that the trial propagator 
must satisfy in order to be the exact one. The method 
would be equivalent to Dyson-Schwinger equations. Of 
course, truncation spoils it and the approximation de¬ 
pends on the accuracy of the polarization function that 
can be evaluated up to second order in any gauge by the 
explicit expressions of the present paper. 

In Landau gauge, Ea. dUUl) can be solved numerically as 
an integral equation for the trial propagator T(p), yeld- 
ing a self-consistent second-order propagator that satis¬ 
fies T 2 (p) = T(p ) because of Ea. (l88l) . The propagator 
satisfies a perfect scaling and can be fitted with good 
accuracy by the simple expression 


T{pe) 


z 

P 2 e + m 2 


(91) 


where Z and m are real parameters. In fact, that explains 
the perfect scaling as any form like that, with just two 
free parameters, can be scaled on top of each other by a 
change of units. The result is very close to that found 
by Feynman gauge in Ref. (4ll| . As shown in Fig.8, after 
renormalization the self-consistent propagator can be put 
on the lattice data, but the agreement is worse than found 
by the method of stationary variance in Fig. 4 and 5 of 
the previous section. Again, the result seems to be almost 
gauge invariant. 

In summary, while the method of stationary variance 
seems to be more reliable than other variational ap¬ 
proaches, other attempts can be devised by the knowl¬ 
edge of the explicit expressions for the polarization up 
to second order. Once optimized, the perturbation the¬ 
ory does not show divergences in the infrared, while the 
ultraviolet ones can be cured by standard regularization 
techniques. The explicit expressions of the present pa¬ 
per hold for any regularization scheme and any choice 
of the gauge parameter. That would suggest a further 
way to optimize the expansion, with the gauge parame¬ 
ter and the renormalization scheme that can be regarded 
as variational parameters [HS3- A comparison between 
the present Landau gauge calculation and previous re¬ 
sults in Feynman gauge [41] shows that the gluon prop¬ 
agator is not too much sensitive to gauge changes and 
that the optimized expansion seems to be approximately 
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gauge invariant after renormalization, which is a desir¬ 
able property from the physical point of view, in qual¬ 
itative agreement with some recent results by Dyson- 
Schwinger equations [H, HH . 


Four-dimensional integrals of functions of the two vari¬ 
ables {k e ■ Pe) and k E are reduced to two-dimensional 
integrals according to 
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Appendix A: Details on the numerical integration 

For d = 4, all numerical integrations have been calcu¬ 
lated as successive one-dimensional integrations by the 
standard Simpson method in the Euclidean space and 
with an energy cutoff p 2 E < A 2 . Four-dimensional in¬ 
tegrals of simple functions of k 2 E are reduced to simple 
one-dimensional integrals before numerical integration, 
according to 


d 4 fc_E 

(2^ 


A{{k E ■ p E ),k%\ = 

y 2 dy 


A „.2 An, rV A2 ~y 2 


4 n 3 


~V A2 ~i 


_A[(xp E ), (a : 2 + y 2 )]dx. 

(A2) 


d A k 


E 


(2tt) 


4 A(k%) — 


1 

87T 2 


A{k 2 )k 3 dk. 


(Al) 


The eight-dimensional integral of the two-loop sunset 
graph (2c) Eg. (17^1) can be written as a four-dimensional 
integral by exact integration of some variables. We no¬ 
tice that each single term contributing in Ea. dT^l) can be 
written as 


J ( 27 t) a J (pe^e^eiIe ■ (k E + pe),Pe ■ k E ) gi (pE,k%,q E ,q E ■ k E ,q E ■ Pe,Pe ■ k E ) (A3) 

where the function /) has one argument less than the function gi , and p is the external momentum. Let us introduce 
the vector V = k E + p E and split the four-vector q as the sum of two orthogonal two-dimensional vectors ( 91 , 92 ) 
and (q x ,q y ) that are orthogonal and parallel to the k — p plane, respectively. Moreover we take the direction q y to 
be parallel to the direction of V. Omitting the variables p 2 E ,k 2 E ,V 2 that are constant in the internal integration, the 
integral in Ea. (|A3l) reads 


f d 4 fc g 1 

/A (2^) 4 2(27t) 4 J 0 


A 2 j-q ^ ny/q 2 -ql 

dq 2 I d(j) / d q y fi{q 2 ,q y ) / _ dq x g t (q 2 , q Xl q y ) 

j-q 


(A4) 


/o 


where tan <f> = 92/91 and q 2 = q E . The angle <j> can be integrated exactly yielding a factor of 27r and denoting by 
gi{q 2 ,q y ) the integrated function 


/V ? 2 -? 2 

9i(q iQy) = I _d q x gi(q ^qx^qy) 


(AS) 


according to Eq. CZD, in Landau gauge we can write the transverse polarization term as 

n f 2 C )(p) = j A $ J o dq2 J d qyKg 2 .gy)J2^ q2 ' q y)9i(g 2 ,g v ) (A6) 


where the function h is 


h{q 2 ,q y ) = T(q E )T{k E )T(k E + q E + p E ) (A7) 

that only depends on q E ■ V ~ q y and 9 2 = q E in the internal integration. The eleven functions gi turn out to be 
polynomials, and the integrated functions <jh can be evaluated exactly so that we are left with a two-dimensional 
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internal integration and the result can only depend on kg and Ue • p- Then, the external integration follows by 
Eq. (IA2l) . yielding a total four-dimensional integration to be evaluated numerically. The eleven terms f i: g iy < 7 ,; follow 

from inspection of Ea. (177ll . Dropping the E in the Euclidean vectors, denoting by s = q 1 2 3 4 5 6 — q 2 and by k x ,k y ,p x ,p y 
the components of Ue, Pe that are parallel to q x and q y , respectively, we can write 


/i=4, 
h = 


1 


h = 
h = 
h = 
h = 
h = 

fs = 

h = 

fw = 
fll = 


q 2 k 2 ’ 
1 

q 2 p 2 ’ 


f + Vq v 


k 2 q 2 {k + q + p) 2 
-ip-k) 


k 2 q 2 p 2 


p 2 (k + p + q) 2 1 

4 <jp-V) 

p 2 (k+p + q) 2 ’ 

2 (P-V) 2 

p 2 (k + p + q) 2 ' 

-2 (p-k) 

p 2 k 2 (k + p + q) 21 
1 

p 2 k 2 q 2 (k + p + q) 2 ’ 
-{k-p){q 2 + Vq y ) 


p 2 k 2 q 2 {k+p + q) 21 
where the functions I n are defined as 

-/l(/u) (2 S^Qyky 

h{k,p) = ( 2 s) ~ 2 
2 


gi = 2s 

32 = h (k,k) 

33 = h(p,p) 

34 = h(k, k) + (V ■ k)I\{k) 

35 = h (p, k ) 

36 = h(p,p) 

37 = h{p) 

3s = 2s 

39 = h(k,p) + (k ■ V)h(p) + (p • V)h(k) + (2s)(p • V)(k ■ V ) 

310 = {p- Vfh{k,k) + h(k,p) + 2 (p- V)I 3 (k,k,p ) 

311 = h(k,p) + (p • V)h(k) 


(AS) 


q^kyPy + -s 2 k x p x 


h{k, k,p) — g s q y (k x p y + 2 k x k y p x ) + (2s)q y (k y p y ) 

Ii(k,p) = (k x p x ) + ^s q y {k x p y + k y p x + 4 k x k y p x p y ) + ( 2 s)q y (k y p y ) 

I- 


(A9) 
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